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Abstract 

3y Monte Carlo simulations in the canonical ensemble, we have studied the magnetic anisotropy in Fe/Dy amorphous multilayers. 
This work has been motivated by experimental results which show a clear correlation between the magnetic perpendicular anisotropy 
and the substrate temperature during elaboration of the samples. Our aim is to relate macroscopic magnetic properties of the 
multilayers to their structure, more precisely their concentration profile. Our model is based on concentration dependent exchange 
interactions and spin values, on random magnetic anisotropy and on the existence of locally ordered clusters that leads to a 
perpendicular magnetisation. Our results evidence that a compensation point occurs in the case of an abrupt concentration profile. 
Moreover, an increase of the noncollinearity of the atomic moments has been evidenced when the Dy anisotropy constant value 
grows. We have also shown the existence of inhomogeneous magnetisation profiles along the samples which are related to the 
concentration profiles. 

Ke^y words: Monte Carlo simulation, Heisenberg model, Ferrimagnetic multilayers. Perpendicular magnetic anisotropy 
PACS: 75.40.Mg, 75.30.Gw, 75.25.+Z, 75.50.Gg 
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^1. Introduction 

Amorphous multilayers of transition metal - rare earth 
(TM-RE) compounds display particular interesting prop- 
erties such as a strong magnetic anisotropy which, in some 
conditions, may be perpendicular to the layers [1 2j. Up to 
now, the accurate origin of this anisotropy is not clearly 
understood and several models have been proposed. These 
models are based on the anisotropic distribution of the 
TM-RE pairs [l], dipolar interactions [3], local structural 
anisotropy |4|5j or single- ion anisotropy ^. However, no 
inodel is able to explain all the experimental features. 

This study has been performed by taking into account re- 
cent experimental results of the structural and of the mag- 
netic properties of amorphous Fe/Dy multilayers |7|8|9j . 
In particular, hysteresis loops measured with the magnetic 
field applied successively in the film plane and perpendicu- 
lar to it show that the perpendicular magnetic anisotropy 
is stronger for the multilayers deposited at a higher tem- 
perature (600K compared to 320K) ^. In this tempera- 
ture range, the structural investigations have shown that 
the multilayer is alternatively composed of layers with a 
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majority of amorphous Fe and of Fe-Dy alloyed layers (Fig- 
ure [1]) [8J. For Ts = 320K, the multilayer is made up of 
roughly pure amorphous Fe and Dy layers with thin inter- 
faces. Above TOOK, the layered structure is destroyed, re- 
placed by mixed grains of a-Fe and Fe2-Dy which leads to a 
planar anisotropy. Altogether, these results show that the 
perpendicular magnetic anisotropy seems to be correlated 
in these samples with the existence of Fe-Dy amorphous 
alloy in the multilayers. 

In this work, we have performed Monte Carlo (MC) sim- 
ulations on amorphous Fe/Dy multilayers in order to in- 
vestigate the influence of the concentration profile on the 
magnetic properties. We have considered a model includ- 
ing perpendicular magnetic anisotropy (on average) on a 
small fraction of Dy sites, the other Dy sites exhibit ran- 
dom anisotropy. Two different concentration profiles have 
been studied: an abrupt profile that is the multilayer is 
made up of pure Fe layers and pure Dy layers (close to the 
experimental profile for Ts = 320K) and an experimental 
profile obtained from tomographic atom probe analysis on 
a multilayer deposited at 570K (Figure [1]). Our goal is to 
investigate by means of MC simulations the influence of the 
concentration profile on the global magnetisation, on the 
global susceptibility, on the sublattice magnetisations and 
on the magnetisation per plane. Since the strength of the 
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Joy-Dy/te = 6.5 (K). 
These concentration rules are valid only when Xpe > 0.4. 
^Fe-Dy IS strongly negative; it is thus responsible for the fer- 
rimagnetic order (with antiparallel Fe and Dy moments) at 
low temperature which has been experimentally observed 
by polarized neutrons reflectivity measurements [9j. We 
have to mention that the Fe spin value also depends on the 
local concentration [lOj: 

SFeiXFe) = LI - 1.125(1 - Xpe) (Xpe > 0.4). 

The other numerical values of the model are those of the 
free ions: 



^Fe 



2; ^Dy = 4/3; ^Dy = 2.5. 



Fig. 1. Concentration profile of an (Fe 3nni/ Dy 2nni) multilayer 
deposited at 570K. The data have been obtained by tomographic 
atom probe measurements 0. 

Dy single site anisotropy is not clearly elucidated, we have 
also investigated the influence of the Dy anisotropy con- 
stant magnitude. Moreover, the fraction of Dy sites with 
perpendicular magnetic anisotropy is not known, so we have 
also let vary this parameter in our simulations to determine 
its influence. The model is presented in section [2] and the 
simulation method is described in section[3l The results are 
discussed in section|4]and a conclusion is given in section[5l 

2. Model 

2.1. Magnetic parameters 

Our model of Fe/Dy ferrimagnetic multilayers consists of 
classical Heisenberg spins located at the sites of a face cen- 
tered cubic lattice with nearest-neighbour exchange inter- 
actions. We have chosen this lattice in order to reproduce 
the experimental observed compact structure [10'11'12]. As 
the multilayers are amorphous, the exchange interactions 
are modulated by a gaussian distribution (Figure [2|) which 
models the distribution of the interatomic distances in real 
samples [13 . On the other hand, the Jpe-Fe and Jpe-Dy ex- 
change interaction values have been extracted from those 
determined by Heiman et al. [lOj on Fe-Dy amorphous thin 
films by mean- field calculations and adjusted by MC sim- 
ulations in order to get the pure amorphous Fe Curie tem- 
perature (270K [To]). So, they linearly depend on the local 
concentration Xpe following: 



Fe-Fe 



JFe-Fe(^Fe)/^B = 77 
-^Fe-Dy(-^Fe)/^B = 8 - 



449(1 -Xpe) (in K) (Xpe > 0.4), 
198(1 - Xpe) (in K) (Xpe > 0.4). 

Concerning the Joy-Dy exchange interaction which is 
known to be much smaller than the others, this latter has 
been taken constant in function of the concentration. Its 
value has been estimated by MC simulations in order to 
provide the pure amorphous Dy Curie temperature (llOK 

my- 
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Fig. 2. Gaussian distribution of the exchange interaction ((a) Fe-Fe 
exchange interaction; (b) Dy-Dy exchange interaction). 

The variation of the exchange energy as a function of 
the Fe concentration for the three different bonds is shown 
in Figure [3l We can see that the exchange energy (in ab- 
solute value) of the Fe-Fe and Fe-Dy bonds are maximum 
respectively at Xpe = 0.8 and Xpe = 0.5 and that the con- 
tribution of the Fe-Dy pairs will play a significant role in 
the magnetic ordering. 

2.2. Concentration profile 

The introduction of realistic concentration profiles along 
the growth direction of the multilayer is essential since 
atomic diffusion seems to have a major influence on the 




Fig. 3. Exchange energy as a function of the Fe concentration for 
the three bonds. 
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Fig. 4. Concentration profile of each specie along the perpendicular 
direction for an abrupt (a) and a TAP (b) concentration profile. 

macroscopic magnetic anisotropy. In this study, we have 
thus considered two concentration profiles: the first one, 
called abrupt^ corresponds to a multilayer made up of pure 
Fe layers and pure Dy layers with an abrupt interface (Fig- 
ure [H^a)); the second profile (Figure [H^b)) is directly ob- 
tained from tomographic atom probe analyses of a multi- 
layer (Fe 3nm/Dy 2nm) built up at 570K [8j (Figure [1]). 
This profile is composed of a rich Fe region (FegoDyio) and 
a wide region in which the concentration varies. As it has 
been previously mentioned, this concentration profile leads 
to the maximum of perpendicular magnetic anisotropy ex- 
perimentally observed by Tamion et al. |8|9j . In the follow- 




Fig. 5. Schematic representation of the ordered Fe-Dy clusters. Fe 
atoms are in blue, Dy atoms in green. The arrows show the local 
anisotropy axes on the Dy sites. 



ing, this profile will be called TAP profile (for tomographic 
atom probe profile). The main difference between the two 
profiles is that the pure Dy region in the multilayer with an 
abrupt profile is replaced by an alloyed region with variable 
local concentration in the case of the TAP profile. 



2.3. Magnetic anisotropy 

In order to investigate the perpendicular magnetic 
anisotropy in the Fe/Dy multilayers, we have focused our 
attention on a local structural anisotropy model [4J which 
seems to be well adapted to understand the referenced 
experimental results. 

This model is based on the central hypothesis that dur- 
ing layer-by-layer growth small Fe-Dy crystallized zones at 
a scale of a few interatomic distances are formed defining on 
average a preferential axis perpendicular to the film plane. 
Indeed, they probably exhibit a structure that is reminis- 
cent of those of the defined TM-RE compounds, whose 
hexagonal symmetry leads to anisotropy axes on average 
perpendicular to the layers. These locally ordered clusters 
consist in our model of 13 atoms (one central Fe atom and 
its 12 nearest neighbours) (Figure [5|). Among the 8 neigh- 
bours which are not in the plane of the central atom, be- 
tween 2 and 4 Dy atoms are randomly distributed in order 
to obtain a cluster concentration close to that of defined 
compounds. Since RE atoms are characterized by a strong 
intrinsic anisotropy due to the particular shape of their 4f 
orbit als, we have thus considered in the Hamiltonian a sin- 
gle site anisotropy term for all Dy atoms. The anisotropy 
axes on the Dy atoms belonging to the clusters are parallel 
to the Fe-Dy bonds (Fe here is the central atom of the clus- 
ter); so the unit vectors z^ of these 4 axes are (±4=, 0,4=) 

or (0, ±4=, 4^). This leads on average to an anisotropy di- 
rection perpendicular to the plane of the layers. For the 
other Dy sites which do not belong to the clusters we have 
considered random anisotropy because of the amorphous 
structure of the multilayers. The energy of the system is 
thus: 



<i,j> i^clusters 

iGclusters 



(1) 



where S^ is the spin of the site z, Jij is the exchange inter- 
action between the nearest-neighbour sites i and j, I^Dy is 
the anisotropy constant on the Dy sites, n^ is an unit vec- 
tor indicating the random anisotropy direction on the Dy 
site i and z^ is an unit vector along the local anisotropy 
axis of a Dy site belonging to the clusters. We also define 
the magnetic susceptibility: 

.2 



X(T) 



{m^)^ - (m) 



^ xTV, 



where (m)^ is the magnetisation defined as 



Vt ~ AT 



N 






1/2. 



As the anisotropy coefficient I^Dy in amorphous multi- 
layers is not accurately determined in the literature, we will 
consider it as a free parameter with the same value for all 
the Dy sites (in the clusters and in the amorphous matrix). 
In the same way, the cluster concentration, which is defined 
as the number of atoms included in the clusters divided by 
the total number of atoms, is a free parameter in the sim- 
ulations since it cannot be evaluated experimentally. But 
in this work, most of the results have been obtained with a 
quite small cluster concentration of 5%. 

3. Numerical method 



exp(— AE'/ZcbT) in order to take into account the ther- 
mal fiuctuations. 

One MC step consists in examining all spins of the sys- 
tem once. At each temperature, 5000 MC steps were per- 
formed to reach the thermodynamic equilibrium, and after- 
wards the physical quantities were measured by averaging 
over the next 10000 MC steps. As we are interested in dis- 
ordered systems, it is necessary to perform simulations on 
several disorder configurations (typically between 20 and 
100) in order to estimate the thermodynamic quantities of a 
macroscopic system. It is important to note that for saving 
CPU time, we have limited our simulations to a double bi- 
layer. In order to match up to the experimental system (Fe 
3nm/Dy 2nm), we have considered 12 Fe and 6 Dy atomic 
planes in each bilayer because the ratio of the atomic Fe 
and Dy radii is around 1.5. Each plane contains 800 atoms 
(which gives a total number of atoms equal to 28800). Since 
we are not interested in the critical behaviour, finite size 
effects have no significant infiuence on the physical proper- 
ties under consideration. 

4. Numerical results 

We present here the magnetic properties of Fe/Dy mul- 
tilayers with the concentration profiles defined previously 
(called abrupt and TAP profiles, see Figure [4j). As a first 
step, since we are interested in the infiuence of the concen- 
tration profile on the magnetic ordering, we have consid- 
ered a multilayer without any anisotropy term. The next 
step of this study is then focused on the effect of the mag- 
netic anisotropy on the Dy atoms in the clusters and in the 
amorphous matrix. 



The numerical method is the importance sampling MC 
method in the canonical ensemble |15|16|17] . The thermo- 
dynamic equilibrium at each temperature, corresponding 
to the minimization of the free energy of the system, is ob- 
tained by the Metropolis algorithm [18] . In this algorithm, 
the spins are examined individually. A site i is randomly 
chosen and a unit vector u^ defined by the random choice 
with uniform distribution of its z-component Ui G [—1,1] 
and an its azimuthal angle (fi G [0, 27r[ is determined. The 
new spin S^ is then given by: 

(^f) = Si^ I- uj cos ifi, 
(^f)'= 5,^1^ sin (^„ 

where Si — Spe or S^y It has to be noted that this spin trial 
rotation procedure is isotropic. Then the energy variation 
A^ associated to this rotation is calculated. The next step 
is the following: 

- if AE < 0, the rotation is accepted; 

- if AE > 0, the rotation may be accepted with a prob- 
ability that is proportional to the Boltzmann factor 



4.1. Influence of the concentration profile without 
anisotropy f^Dy = 0) 

No anisotropy term has been introduced into the Hamil- 
tonian, this latter is then invariant under any rotation of 
all the spins. In these samples, the ground states consist 
of a ferromagnetic Fe layer and a ferromagnetic Dy layer 
which are antiparallel. These ground states, so called fer- 
rimagnetic, will be considered as a reference in the follow- 
ing. The thermal variation of the magnetic susceptibility 
for the whole sample and for each sublattice is shown in 
Figure [6] with the abrupt and TAP concentration profiles. 
We observe that the abrupt profile is characterized by two 
separate magnetic ordering at respectively 285K and 150K 
corresponding to the Fe layer and the Dy layer. The Dy sus- 
ceptibility curve (Figure injb)) exhibits a shoulder at T ~ 
280K related to the polarisation, via the exchange interac- 
tion, of the Dy atoms located in the neighbourhood of the 
interface. For the TAP profile, the two sublattices order at 
the same temperature (325K). The shoulder of the Dy sus- 
ceptibility curve below the peak temperature corresponds 
to the more concentrated Dy planes. The thermal variation 
of the magnetisation of the whole Fe/Dy multilayer (Fig- 
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Fig. 6. Thermal variation of the magnetic susceptibihty for a mul- 
tilayer (Fe 3nm/Dy 2nm) with the 2 studied concentration profiles 
((a) for the whole sample, (b) per sublattice). The curves for the 
TAP profile result from the average over 4 disorder configurations. 

urellKa)) is strongly influenced by the concentration profile 
even if the concentration of each specie is the same in the 
two cases (Xpe = 2/3). 

We note in the case of the abrupt profile the presence of 
a compensation point at T = 175K when the two sublat- 
tice magnetisations are antiparallel with the same modu- 
lus. Concerning the TAP profile, the magnetic ordering is 
continuous and smooth because of the mixing of the two 
species. 

4.2. Influence of the magnetic anisotropy: sperimagnetic 
order 

4.2.1. Metastable states 

The ground states result from the competition between 
exchange interactions and magnetic anisotropy. The ex- 
change interactions lead to a collinear order in these sys- 
tems since there is no frustration whereas the magnetic 
anisotropy favours a angular distribution of the magnetic 
moments. Then the ground states are expected to be sper- 
imagnetic with a cone angle that should increase with the 
anisotropy constant l^Dy 

As it has been described in the section 12.31 due to the 
locally ordered clusters, the average orientation of the mo- 
ments should be perpendicular to the layers. Because of 
the competition between exchange and anisotropy ener- 




Fig. 7. Thermal variation of the reduced magnetisation for a multi- 
layer (Fe 3nm/Dy 2nm) with the 2 concentration profiles ((a) for the 
whole sample, (b) per sublattice). Mferrimagnetic {^) corresponds to 
the ground states with antiparallel Fe and Dy layer magnetisations 
whereas Mgat (b) is the saturated magnetisation of each ferromag- 
netic sublattice. 



gies, numerical convergence towards metastable states at 
low temperature can occur during the simulation, mainly 
in the case of small cluster concentrations and small values 

of L^Dy. 

Figure[8]shows the histograms over 100 chemical disorder 
configurations of the perpendicular magnetisation compo- 
nent of the Dy sublattice for a multilayer with the TAP 
profile at T = IK. The concentration of ordered clusters is 
of 10% and I^Dy/fe = lOK. Without any random magnetic 
anisotropy (RMA) on the Dy atoms of the matrix (Figure 
[51^ a)), we observe a sharp distribution at |Mgy/MDy| ~ 1 
which means that the magnetisation of the Dy sublattice is 
perpendicular to the layers for all disorder configurations. 
The introduction of RMA on the Dy atoms of the ma- 
trix leads to a broader distribution P(|M^ /MDy|) (Fig- 
ure [8jb)). In particular, some disorder configurations ex- 
hibit a magnetisation which is clearly non perpendicular to 
the layers. This phenomenon is less pronounced when the 
anisotropy constant I^Dy increases (Figure [8l^c)). Anyway, 
this problem has been overcome, by decreasing the cooling 
rate during the simulation and also increasing the number 
of MC steps at each temperature. 
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Fig. 8. Histograms over 100 disorder configurations at T = IK of 
the perpendicular magnetisation component of the Dy sublattice 
measured with the TAP profile, a concentration of clusters of 10% 
((a) without any RMA and Dj^y/kB = lOK, (b) with RMA on the 
matrix Dy atoms and I^Dy/^B = lOK, (c) with RMA on the matrix 
Dy atoms and DDy/^B = 50K). 
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Fig. 9. Thermal variation of the reduced magnetisation for a mul- 
tilayer (Fe 3nm/Dy 2nm) with an abrupt concentration profile con- 
taining 5% of clusters for different values of the Dy anisotropy con- 
stant -Doy ((a) total magnetisation, (b) Fe and Dy sublattice mag- 
netisations). 

4.2.2. Influence of the concentration profile for different 
I^Dy values 

Here we have studied the influence of the magnetic 
anisotropy constant on the Dy sites in a multilayer contain- 
ing 5% of locally ordered clusters which means that only 
3.5% of Dy atoms exhibit uniaxial anisotropy (on average). 

• Abrupt profile 
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Fig. 10. Schematic representation of the magnetic order at low tem- 
perature in a multilayer with an abrupt concentration profile for 3 
different anisotropy constant values I^Dy 



We firstly presente the results concerning the abrupt con- 
centration profile. Figure [9] shows the thermal variation 
of the global magnetisation and of the Fe and Dy sublat- 
tices for different values of I^Dy • We observe a strong influ- 
ence of I^Dy on the low temperature magnetisation (T < 
^comp = 175K) (Figure [91(a)), especially on the Dy sublat- 
tice, its influence on the Fe sublattice being imperceptible 
(Figure [91(b)). This phenomenon can be understood by the 
fact that the magnetic ordering process above Tcomp is es- 
sentially due to the Fe magnetic moments which do not 
display any magnetic anisotropy. Below Tcomp, both the 
low-temperature global magnetisation and the Dy sublat- 
tice magnetisation decrease when the anisotropy constant 
increases. This is due to the increase of the non collinearity 
of the Dy moments, that is a sperimagnetic order already 
described in the literature in similar systems |19|2Qj . These 
magnetic configurations are qualitatively shown in Figure 
[To] for three values of D^y . The angular distribution of Dy 
moment is broader in the core of the Dy layer whereas it is 
the contrary for the Fe moments which are more misaligned 
close to the interface because of Fe-Dy coupling. 

• TAP profile 
The influence of the concentration profile on the sperimag- 
netic order has been investigated, using the previously de- 
fined TAP profile. Figure [TTl shows the thermal variation of 
the global magnetisation and of the Fe and Dy sublattices 
for a concentration of ordered clusters of 5% and different 
values of I^Dy • 

As in the case of the abrupt profile, we observe a decrease 
of the global magnetisation of the multilayer when Dj^^y 
increases but this decrease is much less marked than in the 
previous case. Indeed, we note that the magnetisation of 
the Dy sublattice is less influenced by I^Dy whereas a small 
effect on the Fe sublattice can be seen when I^Dy/^B = 50K. 
It means that as the layers are more mixed, the influence 
of the RMA on the Dy atoms is partially transmitted to 
the Fe sublattice which displays also a low temperature 
non-collinear magnetic order (Figure [T2|). We can see that. 
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Fig. 11. Thermal variation of the reduced magnetisation for a multi- 
layer (Fe 3nm/Dy 2nm) with a TAP concentration profile containing 
5% of clusters for different values of the Dy anisotropy constant I^Dy 
((a) total magnetisation, (b) Fe and Dy sublattice magnetisations). 
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Fig. 12. Schematic representation of the magnetic order at low tem- 
perature in a multilayer with a TAP concentration profile for 3 dif- 
ferent anisotropy constant values I^Dy 

at least for I^Dy/^B = 50K, the angular distribution of 
the Fe moments broadens while those of the Dy moments 
sharpens. 

4.2.3. Influence of the concentration profile on the 
magnetisation profiles (D^y/kB = bOK) 



In order to get a quantitative analysis of the sperimag- 
netism, we have calculated the magnetisation for each plane 
and the histogram of the perpendicular component of the 
magnetic moments for each plane and each sublattice. Since 
the sperimagnetism is more pronounced for I^Dy/^B = 
50K, we have chosen this case. The results obtained at T = 
IK are shown in Figures [T3l and [T4l for the abrupt and TAP 
profiles respectively. 

For the abrupt profile, we clearly see the evolution of the 
magnetisation per plane which is not homogeneous inside 
the multilayer (Figure [T3|) . The decrease of the magnetisa- 
tion is related to the broadening of the distribution of the 
perpendicular component of the atomic moments as men- 
tioned previously. More precisely, far from the interface, the 
distribution is very broad in the Dy atomic planes (Figure 
ITSTa)) due to the large value of the ratio l^Dy/-^Dy-Dy = 
7.7. At the interface, the distribution width is minimum in 
the Dy layer (Figure [T3lf b)) whereas it is maximum in the 
Fe layer (Figure [TS^c)) because of the strong Fe-Dy cou- 
pling. Finally, the distribution is close to a Dirac peak in the 
core of the Fe layer since there is no magnetic anisotropy 
(Figure [md)). 

In the case of the TAP profile (Figure [T4j), the variation 
of the magnetisation per plane for each sublattice is dif- 
ferent from the previous case since it depends on the local 
concentration. We observe that the magnetisation of the 
Fe sublattice reaches a maximum for the high-concentrated 
Fe planes (Xpe = 0.9) whereas for the Dy sublattice, the 
magnetisation per plane exhibits a maximum for a plane 
with Xpe close to 0.4. This agrees with the distribution of 
the perpendicular component of the magnetic moments for 
the Dy sublattice which is significantly sharper in planes 
with Xpe = 0.4. Contrarily, for the Fe sublattice, the distri- 
bution is sharper in the rich Fe planes. In order to explain 
these features, we have reported in Table [H the magnitude 
of each energy term (per atom) defined as: 
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where X^^^^^^ is the Dy atomic fraction in the clusters. In 
the 24*^ plane (rich-Fe plane), the coupling between the 
Dy and Fe sublattices is weak, then the Fe moments are 
almost collinear whereas the Dy moments are strongly out 
of line due to the large ratio £^R^A/^e^xch ^^ + ^e^xc"h^^ ^ 
5.2. Concerning the 16*^ and 13*^ atomic planes, the Dy- 
Fe sublattice coupling is strong and the possible speri- 
magnetism of the Dy sublattice (due to the very large 
^RMA value) could be transmitted to the Fe sublattice. 
The main difference between these two planes is the ratio 



E, 



Dy 
RMA 



/(^exch 



Dy-Dy 



^exdi ^) which is about 1.9 and 2.1 in 
the 13*^ and 16*^ atomic planes respectively. This explains 
why the magnetic configuration is more collinear in the 13*^ 
plane than in the 16*^ plane. The comparison of these two 
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Fig. 13. Reduced magnetisation of each atomic plane of a multilayer at low temperature (T = IK) with the abrupt concentration profile 
containing 5% of clusters and i^Dy/^B = 50K. Distribution of the perpendicular magnetic component in the Dy layer and in the Fe layer 
((a) Dy - core plane; (b) Dy - interface plane; (c) Fe - interface plane; (d) Fe - core plane). 
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Fig. 14. Reduced magnetisation of each atomic plane of a multilayer at low temperature (T = IK) with the TAP concentration profile 
containing 5% of clusters and I^Dy/^B = 50K. Distribution of the perpendicular magnetic moment component in the Dy sublattice and in 
the Fe sublattice ((a) plane 16 - X^e = 0.50; (b) plane 24 - Xpe = 0.90) 



values with those of the 24*^ plane accounts for the more 
marked sperimagnetism of the Dy sublattice in this latter. 

Table 1 

Random magnetic energy and exchange interaction energy for 3 

typical atomic planes for the TAP profile. 



plane Xpe 


<^A (K) 


^Dy-Dy 
exch 


(K) 


^r^^^ (K) 


<:cr (K) 


13 0.65 


102.34 


4.97 




49.25 


49.34 


16 0.50 


146.88 


10.16 




61.14 


21.77 


24 0.90 


29.30 


0.41 




5.24 


86.29 



These results evidence that the magnetisation profile 
along the multilayer is not homogeneous, which has already 
been observed experimentally on these samples from po- 
larized neutrons reflectivity measurements [9j. We would 
like to emphasize that unlike what is mentionned in sev- 
eral papers [21 22 23j, our results indicate that the atomic 
moment distribution around the mean magnetisation di- 
rection in each plane is not uniform but rather gaussian for 
the two concentration profiles. 

4.2.4. Influence of the cluster concentration 

The influence of the cluster concentration also has been 
investigated for I^Dy/te = 50K in the case of the TAP 
profile. 
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Fig. 15. Thermal variation of the reduced magnetisation for a mul- 
tilayer (Fe 3nm/Dy 2nm) with a TAP concentration profile for 
D-Dy/k-Q = 50K for different cluster concentrations ((a) total and 
sublattice magnetisation, (b) magnification of the sublattice mag- 
netisation). 

Figure [15] shows the thermal variation of the reduced 
magnetisation for such a multilayer for different cluster con- 



centrations. The increase of the cluster concentration from 
5% to 20% has no influence on the global magnetisation: 
the three curves for each concentration are superimposed 
at each temperature (Figure fTSlf a)). Looking more precisely 
at the sublattice magnetisations at low temperature (Fig- 
ure [T5jb)), we note a very small increase of each sublat- 
tice magnetisation as the cluster concentration rises. Quite 
surprisingly, the cluster concentration has no significant ef- 
fect on the angular distribution of the moments unlike the 
anisotropy constant to which the magnetic order is very 
sensitive. This result can be interpreted from the strong 
influence of the exchange interactions in these multilayers 
which induce a global behaviour of the magnetic moments 
as previously seen. 

5. Conclusion 

In order to qualitatively describe the macroscopic per- 
pendicular anisotropy, we have used in our model a local 
structural anisotropy model which induces a perpendicular 
magnetisation even for very small cluster concentrations. 
This study has shown using a numerical investigation the 
influence of the concentration profile on the thermal varia- 
tion of the magnetisation and on the magnetisation profile 
for several values of I^Dy For the abrupt profile, a compen- 
sation point can be observed and sperimagnetism occurs 
essentially in the Dy layer. Concerning the TAP profile, the 
sperimagnetism on the Dy sublattice is not homogeneous 
and is less pronounced in the rich Dy plane with Xoy ~ 
0.60. In a near future, we plan to perform simulations of 
hysteresis loops in order to explain experimental data by 
suggesting rotation mecanisms for the magnetic moments 
in function of the concentration profile, the anisotropy con- 
stant I^Dy and the cluster concentration. 
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